# Additional R code for lecture 5

# Fig. 3
apply(fac.vals2[1:6,],1,function(x) paste(x[1:2],sep='',collapse=''))->mylabs
oldmar<-par("mar")
par(mar=c(2.1,6.1,4.1,2.1))
plot(c(0.7,5.5),c(0.7,5.5),type='n',axes=F,xlab='',ylab='')
 axis(3,at=1:5,mylabs[2:6])
  axis(2,at=5:1,mylabs[1:5], las=2)
  box()
  k<-1
  for(j in 5:1){
    for(i in (6-j):5){
         points(i,j,pch=16,cex=.7,col=mycol[(out.mat[k]<.8)+1])
         text(i,j,round(out.mat[k],3),pos=3,col=mycol[(out.mat[k]<.8)+1])
         k<-k+1}}
 mtext(expression(bold('Mean 1')),side=2,line=4, cex=1.1)
 mtext(expression(bold('Mean 2')),side=3,line=2.5, cex=1.1)

# Fig. 4b
plot(c(0,3),c(.5,4.5),type='n',axes=F,xlab='',ylab='')
 k<-1
 for(j in 1:4){
 for(i in 1:2){
 points(i,j,pch=16,cex=3,col=k)
 text(i,j,k,pos=3,cex=1.2,offset=0.9)
 k<-k+1}}

# Fig. 4c
plot(c(0,3),c(.5,4),type='n',axes=F,xlab='',ylab='')
 loc<-seq(1,3,.5)
 for (i in 1:6) {
 segments(1,loc[i],2,loc[i],lty=i,lwd=2)
 text(1,loc[i],i,pos=2,cex=1.2)
}